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ABSTRACT 

We report additional photometric CCD observations of KPD 0422+5421, a binary 
with an orbital period of 2.16 hours which contains a subdwarf B star (sdB) and a 
white dwarf. There are two main results of this work. First, the light curve of KPD 
0422+5421 contains two distinct periodic signals, the 2.16 hour ellipsoidal modula- 
tion discovered by Koen, Orosz, & Wade (1998) and an additional modulation at 7.8 
hours. This 7.8 hour modulation is clearly not sinusoidal: the rise time is about 0.25 
in phase, whereas the decay time is 0.75 in phase. Its amplitude is roughly half of 
the amplitude of the ellipsoidal modulation. Second, after the 7.8 hour modulation is 
removed, the light curve folded on the orbital period clearly shows the signature of 
the transit of the white dwarf across the face of the sdB star and the signature of the 
occultation of the white dwarf by the sdB star. We used the Wilson-Devinney code to 
model the light curve to obtain the inclination, the mass ratio, and the f2 potentials, 
and a Monte Carlo code to compute confidence limits on interesting system parame- 
ters. We find component masses of M sdB = 0.36^°;^ M and M WD = 0.47+°;}g M 
(Mtotai = O.86±°;35M , 68 per cent confidence limits). If we impose an additional 
constraint and require the computed mass and radius of the white dwarf to be con- 
sistent with a theoretical mass-radius relation, we find M s( jb = 0-51lio oso -^0 an< ^ 
Jl%D = 0.526lg Mq (68 per cent confidence limits). In this case the total mass of 
the system is less than 1.4 M Q at the 99.99 per cent confidence level. We briefly dis- 
cuss possible interpretations of the 7.8 hour modulation and the importance of KPD 
0422+5421 as a member of a rare class of evolved binaries. 

Key words: binaries: close - stars: individual: KPD 0422+5421 - stars: white dwarfs 
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1 INTRODUCTION 

KPD 0422+5421 was first identified in the Kitt Peak Downes 
survey of blue stars at mid to low galactic latitudes (Downes 
1986). It was classified as a hot, hydrogen- rich subdwarf B 
(sdB) star. The binary nature of this star was discovered 
by Koen, Orosz, & Wade (1998, hereafter KOW) who de- 
termined an orbital period of P = 0.0901795 ± (3 x 10 -7 ) 
days (2.16 hours) using spectroscopic and photometric ob- 
servations. KOW showed that the companion star of the 
sdB is a white dwarf, and determined masses of M s( jb = 
0.72 ± 0.26 Mq and Mwd = 0.62 ± 0.18 M Q for the sdB star 
and white dwarf, respectively. KPD 0422+5421 has one of 
the shortest known orbital periods for a detached binary, 
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and is one of only a few known binaries which contain an 
sdB star paired with a white dwarf. This binary has had 
an interesting past, given that it contains two evolved stars 
in a close orbit, and will have an interesting future, given 
the relatively short time-scale for merger due to the loss of 
orbital angular momentum via gravitational wave radiation 
(~ 1.5 x 10 8 yr), which is comparable to the core He burning 
lifetime of an sdB star (Dorman, Rood, & O'Connell 1993). 

The U and B light curves of KPD 0422+5421 presented 
in KOW are ellipsoidal, with amplitudes of ~ 0.02 mag. 
Ellipsoidal model fits to the light curves using the Wilson- 
Devinney (1971, hereafter W-D) code indicated a high in- 
clination (78°) which should produce an eclipsing geometry. 
The signature of the white dwarf transit across the face of 
the sdB star was predicted to be ~ 0.006 mag deep, last- 
ing ~ 27° in phase (w 10 minutes of time). The presence 
of eclipses would be important since one could in principle 
then determine the geometry of the system reasonably well. 
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However, the U and B light curves of KOW, which were ob- 
tained using high-speed photomultipliers, were sufficiently 
noisy that this subtle signature could not be clearly seen. 
We therefore obtained additional photometric observations 
using a CCD with the goal of obtaining a precise eclipse 
profile. We discuss below our observations, data reduction 
and analysis, and our results. 

2 OBSERVATIONS AND DATA REDUCTION 

We observed KPD 0422+5421 December 27-30, 1998 (all 
dates are UT) using the Kitt Peak 2.1 metre telescope and 
the T1KA 1024 x 1024 CCD. The CCD was binned 2x2 
on chip, yielding a readout time of « 11 seconds and a 
scale of 0'.'6 per pixel. Nearly all of the observations were 
made with a standard Kron R filter (we used the KPNO 
"Harris set" UBVRI filters). The exposure times ranged 
from 20 to 40 seconds. Care was taken to place the im- 
age of KPD 0422+5421 on the same region of the detector 
in order to minimize possible systematic errors introduced 
by flat-fielding errors. The observations on the first night 
were mostly taken through thin cirrus, while the other three 
nights were clear. The seeing was typically 1.0 to l'. 2, al- 
though on some occasions the seeing was as low as 0'.'75. We 
obtained a total of 2676 i?-band images of KPD 0422+5421 
covering 9.17 hours, 8.83 hours, 8.75 hours, and 9.83 hours 
on December 27 through December 30, respectively. On De- 
cember 28 we also obtained images of a total of 38 standard 
stars in four different fields (RU 149, RU 152, PC 0918+029, 
PC 1047+003) from the list of Landolt (1992) in the U, B, 
V, R, and / filters (airmass range 1.19 to 1.97), as well as 
three sets of images of KPD 0422+5421 in the same filters 
(airmasses of 1.44, 1.46, and 1.63). 

Standard IRAF^ tasks were used to remove the elec- 
tronic bias and to perform the flat-fielding corrections. The 
programs DAOPHOT He, ALLSTAR and DAOMASTER 
(Stetson 1987; Stetson, Davis, & Crabtree 1991; Stetson 
1992a,b) were used to compute the photometric time se- 
ries of KPD 0422+5421 and all of the field stars. The 
DAOPHOT and ALLSTAR codes were used because KPD 
0422+5421 has a faint neighbour star « 2'.'3 to the south 
(AR — 5.1) and because many of the bright field stars 
also have near neighbours. We used an iterative proce- 
dure to arrive at the final light curve. First, an automated 
DAOPHOT/ ALLSTAR script was used to determine the 
positions and instrumental magnitudes of all the stars for 
each frame. Then using the stellar positions determined 
above, the DAOMASTER code was used to compute the 
mean pixel offsets for each image with respect to the first 
one. The pixel shifts were rounded to the nearest integer (to 
avoid interpolation in the shifting routine) and the images 
were then aligned so that they all had a common coordi- 
nate system. A deep "master" image was made by combin- 
ing 101 images where the seeing was better than « Cf.'9, 
and DAOPHOT and ALLSTAR were used to reduce this 
master image interactively. In the end we selected nine- 
teen relatively bright and isolated field stars which gave 
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the best fits to the empirical point spread function (PSF). 
The final "master photometry list" for this deep image had 
838 stars, compared to the w 300 stars typically found 
by the daofind routine on any given single image. Next, 
the aligned images were reduced using another automated 
DAOPHOT/ ALLSTAR script. In this case, the stars from 
the master photometry list were used in place of the ini- 
tial stellar list generated by daofind for each image and the 
same 19 stars were used to fit the PSFs for each image. We 
found that the script using the larger master photometry list 
and common PSF stars gave much better results than the 
simple script used for the first iteration. Finally, the DAO- 
MASTER code was used to robustly compute the mean 
magnitude offsets for each image and correct the magnitude 
system of each image to the magnitude system of the first 
image by adding a constant. 

The observations of the Landolt standard stars were 
used to calibrate the DAOMASTER instrumental magni- 
tudes to the standard system. We used the IRAF task phot 
to perform aperture photometry on the standard stars using 
a sequence of different sized apertures of radius 3 to 16 pixels 
(l'.'8-9'.'6). The IRAF implementation of the DAOGROW 
algorithm (Stetson 1990) was then used to fit the curve- 
of-growth function for each image and compute the opti- 
mal instrumental magnitudes for the r = 16 pixel aperture. 
The resulting list of instrumental magnitudes for the stan- 
dard stars was used to define these transformations from 
the IRAF instrumental magnitude system to the standard 
system: 

mu = (3.400 ± 0.044) + U + (0.415 ± 0.028)X[/ 

-(0.151 ± 0.077) (U-B) 

+(0.057 ±0M3)(U -B)Xu 
m B = (1.425 ± 0.034) + B + (0.236 ± 0.022)X S 

-(0.023 ± 0.058)(S - V) 

-(0.049 ± 0.039) (B - V)X B 
m v = (1.116 ±0.027) + V + (0.127 ± 0.017)AV 

±(0.045 ±0.046)(B - V) 

-(0.018 ±0.030) (B - V)X V 
m R = (1.065 ± 0.018) +R + (0.053 ± 0.011)X fl 

-(0.004 ± 0.052) (V - R) 

-(0.019 ±0.034)(V - R)X R 
mi = (1.630 ±0.036) +1+ (0.014 ± 0.023)A+ 

— (0.176 ± 0.112)(i? — I) 

±(0.079 ± 0.075)(i?-I)A+ (1) 

where the lower case ra's are the instrumental magnitudes 
and where X is the airmass. Phot was also used in a similar- 
manner on the December 28 KPD 0422±5421 images (the 
three sequences of U, B, V, R, and I exposures). We pre- 
pared these 15 images by using DAOPHOT to remove every 
star except the image of KPD 0422±5421 and eight bright 
field stars. The transformation equations were inverted to 
derive the calibrated magnitudes of KPD 0422±5421 and 
the eight field stars. Finally, the mean difference between 
the DAOMASTER magnitudes and the standard R mag- 
nitudes for the eight stars was computed and this offset was 
applied to every remaining star to arrive at the final R mag- 
nitude light curves. 

On December 28.3934 (UT) we found the following 
colours and magnitudes for KPD 0422±5421: V = 14.66 ± 
0.02, B — V — 0.18 ± 0.02, U - B = -0.61 ± 0.04, 
V - R = 0.11 ± 0.02, and R - I = 0.11 ± 0.03. For com- 
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Figure 1. Top panels: the time series photometry of KPD 0422+5421 (i?-band). In addition to the ellipsoidal modulation at 2.16 hours, 
a signal with a much longer period is evident. Bottom panels: the light curves of two nearby comparison stars, plotted using the same 
range in the y-axis as the top panels. 



parison, Downes (1986) measured V = 14.66, B — V = 0.20, 
and U — B — —0.65 (no errors were given). 



3 THE PHOTOMETRIC TIME SERIES 

Figure [j] shows the light curves of KPD 0422+5421 from each 
night (top) and the light curves of two bright field stars. The 
standard deviations in these light curves are 0.009 mag for 
the brighter comparison star and 0.011 mag for the fainter 
comparison star. The first observations on the fourth night 
were taken in bright twilight, and as a result the photome- 
try from that time is rather noisy. We therefore excluded the 
data from the first hour of the fourth night (our main results 
are not changed significantly when these data are included 
in the analysis). One can plainly see the 2.16 hour ellipsoidal 
modulation in the light curves of KPD 0422+5421. Surpris- 
ingly, there is a much longer period signal evident in the 
light curves. This signal is especially evident on the third 
night. We argue below that this signal is real and not due to 
an artifact of the observations or of the data analysis. What- 
ever its origin, this extra signal cannot be neglected since its 
amplitude is roughly half that of the ellipsoidal component. 
We therefore must remove it before the ellipsoidal light curve 
can be properly analyzed. 



3.1 Decomposition of the light curve 

We used an iterative procedure to decompose the light curve 
into the ellipsoidal component and the longer-term compo- 
nent. The magnitudes were converted to a linear intensity 
scale and normalized so that an R magnitude of 14.55 cor- 
responds to an intensity of 1.0. There are six steps in each 
iteration: (1) The linearized data were folded on the photo- 
metric ephemeris of KOW and binned into 100 phase bins. 
(2) A W-D model was fitted to the folded light curve using 
a Levenberg-Marquardt optimization procedure (the W-D 
parameters determined by KOW were used for the initial 
estimate). (3) The fitted W-D model was "unfolded" using 
the KOW ephemeris to yield a curve giving the intensity as 
a function of time and subtracted from the linearised time 
series data. (4) The residuals were searched for periodicities 
by computing the variance statistic of Stellingwerf (1978) for 
trial periods between 0.01 and 0.5 days. (5) The residuals 
were phased folded on the significant period found in step 
(4), binned into 250 phase bins, and characterized by a fit 
using a 2 piece cubic spline. (6) The spline fit was unfolded 
and the resulting curve was used to detrend the time series 
data, and the process was started again at step (1) using the 
detrended data. 

The results of the period search from step (4) are dis- 
played in Figure |^. There is a local minimum at a trial pe- 
riod of 0.325 days (7.8 hours), meaning there is a coherent 
modulation at that period (the local minima at 0.245 and 
0.481 days are one cycle per day aliases of the 0.325 day pe- 
riod) . The trial period of the local minimum did not deviate 
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Figure 2. The © statistic of Stellingwerf (1978) vs. trial period 
in days for the KPD 0422+5421 light curve with the ellipsoidal 
modulation removed (solid line). There is a local minimum at a 
period of 0.325 days (the minima at 0.245 and 0.481 days are one 
cycle per day aliases). The G statistic for the comparison star 
with R = 14.715 (see Figure 1) is shown as the dash-dotted line. 

from 0.325 days after the first iteration. Figure |^ shows the 
residuals folded on the 7.8 hour period and binned into 250 
equal phase bins (the phase was adjusted to give the mini- 
mum at phase 0), and the fit using a two piece cubic spline 
(solid line). There is a clear non-sinusoidal modulation with 
an amplitude of 1 per cent of the peak flux. We found that 
the character of this 7.8 hour modulation did not noticeably 
change after the first iteration. 

The ellipsoidal light curve arrived at by folding and bin- 
ning the detrended data did not noticeably change after the 
second iteration. Figure ^| shows the final ellipsoidal light 
curve, converted back into the R magnitude scale, and the 
best-fitting W-D model (i = 89°, M wo /M BdB = 1.15, see 
below for details of the fitting procedures). The signature 
of the white dwarf transit across the face of the sdB star is 
clearly seen at phase 0.0. Surprisingly, a second narrow dip 
in the light curve due to the total eclipse of the white dwarf 
by the sdB star can be seen at phase 0.5. The relatively large 
depth of this feature (w 0.005 mag) indicates that the white 
dwarf is much brighter than supposed by KOW: the new W- 
D fit gives a blackbody temperature of T 2 = 24, 300±2600 K, 
compared with the value of T2 ~ 4000 K derived by KOW. 
The properties of the white dwarf are discussed further be- 
low. 

3.2 Possible systematic errors 

There is clearly an extra periodic signal in the light curve 
of KPD 0422+5421. Removing this signal gives a cleaner 
ellipsoidal light curve. Is the 7.8 hour modulation intrinsic 
to KPD 0422+5421 or is it an artifact of the observations or 
reductions? We consider possible sources of systematic error 
below. 

One could have errors in photometry owing to changes 
in seeing and/or changes in the stellar positions on the CCD 
(see Frandsen et al. 1996). However, since the DAOPHOT 
code uses an analytic model combined with an empirical 
look-up table of residuals when it fits the PSF to each im- 
age, we believe it can robustly handle changes in the PSF 
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Phase (P=0.324B, T =2,451, 174.35348) 
Figure 3. The residuals are shown folded on a period of 7.8 hours 
and binned into 150 phase bins, where the points correspond to 
the median value within that bin. The error bars represent the 
error of the mean within each bin. The solid line is a two piece 
cubic spline fit. Two cycles are shown for clarity. Note that the 
y-axis is a linear flux scale. 

due to guiding errors or changes in seeing. Furthermore, the 
ALLSTAR program recomputes the stellar centroids at each 
iteration based on the most recent PSF. Thus we do not ex- 
pect errors in the photometry due to seeing changes or shifts 
in the stellar positions. Also, one would expect to see spuri- 
ous signals in the light curves of other stars if drifts in the 
stellar positions or changes in seeing affected the photome- 
try. However, no spurious signals of comparable amplitude 
are evident in the light curves of any of the bright field stars. 

The differential light curves produced by the DAO- 
MASTER code were made with the assumption that the 
extinction (either due to clouds or to changing airmass) is 
independent of the colour of the star. If the extinction did 
depend on the colour of the star, then the light curve of 
KPD 0422+5421 could show a spurious trend since KPD 
0422+5421 is bluer than its immediate neighbours (~ 0.5 
mag bluer in V — R than the bluest field stars). We do not 
believe our photometry suffers from this problem since the 
coefficients of the "colour- airmass" terms in the transforma- 
tions (equation 1) are all consistent with zero. Several of the 
standard stars we observed are bluer than KPD 0422+5421 
(e.g. RU 149 with B—V = -0.129 and U-B = -0.779, Lan- 
dolt 1992) and were observed to fairly high airmass (« 2), 
so we would be able to measure colour-airmass terms that 
were significantly different from zero. As an experiment, we 
considered only the KPD 0422+5421 observations obtained 
at airmass less than 1.4 (1977 images) and analyzed the 
residual light curve (i.e. the time series data with the el- 
lipsoidal component removed). The residuals still showed a 
periodicity at 7.8 hours and the residuals folded on 7.8 hours 
closely resembled the folded curve shown in Figure ^. Thus 
the 7.8 hour modulation does not seem to be an artifact of 
colour-dependent extinction. 

None of the other bright field stars show any apparent 
coherent periodicity. We show the O statistic for a com- 
parison star in Figure 2. The minimum O value is about 
0.89. Furthermore, the O values are greater than 0.95 at 
trial periods near 7.8 hours. Thus the fact that only KPD 
0422+5421 shows an apparent coherent periodicity is an in- 
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Figure 4. The detrended data (i.e. the original time series with the 7.8 hour modulation shown in Figure ^ removed) are shown folded 
on the orbital period (0.0901795 days) and binned into 100 bins, where the points correspond to the median value within that bin. The 
error bars represent the error of the mean within each bin. The solid line is the best-fitting W-D model (see the text). One can clearly 
see the signature of the transit of the white dwarf in front of the sdB (phase 0.0) and the signature of the occultation of the white dwarf 
by the sdB (phase 0.5). Two cycles are shown for clarity. The y-axis here is a magnitude scale. 



dication that the cause is not instrumental, and would have 
to be related to the observing conditions (for example the 
potential systematic error associated with the blue colour 
of KPD 0422+5421 discussed above). However, if the extra 
signal in the KPD 0422+5421 light curve is indeed an ar- 
tifact of the observing conditions and/or procedure, then it 
is not obvious how this signal could have a periodicity at 
7.8 hours. One might be able to invent a periodicity at 7.98 
hours (one-third of a sidereal day) if the signal were due to 
the rapid change of airmass as the source sets. The fact that 
the periodicity of the signal is not obviously related to the 
rotation of the Earth is further evidence that it is intrinsic 
to KPD 0422+5421 and not an artifact. 



4 ANALYSIS OF THE REFINED 
ELLIPSOIDAL LIGHT CURVE 

4.1 Model outline 

With the 7.8 hour modulation removed from the time se- 
ries data, it was possible to properly analyze the ellipsoidal 
light curve shown in Figure ^| We modelled the folded light 
curve using "mode 2" of the W-D code, which is the usual 
mode used for detached binaries. We assumed a circular or- 
bit and synchronous rotation. The polar temperature of the 
sdB star was fixed at the spectroscopically determined value 
of 25,000 K (KOW), and we used coefficients appropriate for 



the linear cosine limb darkening law taken from the tables 
of Van Hamme (1993). The gravity darkening exponents for 
both stars were set to 1, the usual values for stars with ra- 
diative envelopes (e.g. Wilson & Devinney 1971 and cited 
references) . [In this particular case, the maximum tempera- 
ture change over the star is about 4 per cent. Although we do 
not expect a gravity darkening exponent different than 1, our 
fits are not sensitive to small changes (jS 10 per cent) in this 
parameter.] Similarly, the bolometric albedos for reflective 
heating and re-radiation were also set to 1. (We note that 
there is virtually no "reflection effect" in KPD 0422+5421 
since the sdB star is more than 5.5 mag more luminous than 
the white dwarf. Hence our fits are not sensitive to the ex- 
act details of the reflection effect computations.) The free 
parameters in the model were the inclination i, the mass ra- 
tio Q = MwD/M B dB, the surface potentials Q B dB and SJwd, 
and the polar temperature of the white dwarf T po i e (WD). 

The "grid size integers" used were n\ = 220 and 712 = 
20. A large value of ni was needed to obtain a smooth light 
curve near phase 0.0, when the white dwarf passes in front 
of the sdB star. Thus it was necessary to change an input 
format statement in the W-D code. In that same vein, we 
also altered the appropriate input format statement to allow 
for fiwD values greater than 99.9999. 
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Figure 5. A contour map of the \ 2 values is shown (for 100 data bins). The filled circle marks the point with the lowest \ 2 value 

(Xmin = 94 - 84 at * = 89 -°° and Q = L15 )- The contour levels shown are for X 2 = xLn + Xmin + 2 - 30 > Xmin + 6 - 17 - Xmin + 1L8 > and 
Xmin IS- 4 - The latter four contour levels correspond to confidence levels of 68.3, 95.4, 99.73 and 99.99 per cent for two parameters of 
interest. W-D solutions with mass ratios below the dash-dotted line have M tota i > 1.4 Mq, while solutions with mass ratios below the 
dotted line have M WD > 1.4 Mq. 



Table 1. Parameters for KPD 0422+5421 



Parameter 



Value with no 
added constraints 



Value with white 
dwarf mass-radius constraint 



Reference 



Orbital Period (days) 
T (photo) 

.fCjdB velocity (km s _1 ) 
orbital separation (-Rq) 



0.09017945 ± (1.2 x 10~ 7 ) 

HJD 2,451,174.493408 ± 0.000057 

237 ± 18 

n 81+ - 11 
U - S1 -0.15 



same 
same 



0.853 



+0.002 



This work 
This work 

KOW 
This work 



i (degrees) 

Q = M WD /M sdB 

f!wD 

Tpoie(sdB), spectroscopic (K) 
Tp ie(WD), W-D model (K) 



> 84.4 

c;+0.11 
-0.34 
-+0.27 
-0.82 
3+7.50 
-28.97 

25, 000 ± 1500 
22, 600 ± 2600 



Lis: 

5.75j 
71.73"' 



> 85.6 

7 +0.09 
-0.14 
-+0.24 
-0.30 

62.81±1;|| 



1.07: 
5.55^ 



This work 
This work 
This work 
This work 

KOW 
This work 



M sdB (Mq) 
A/WD (Mq) 
total mass (Mq) 
R s dB, polar (Rq) 
Rwd, polar (Rq) 
log9adB, W-D model 
log9 s dB> spectroscopic (cgs) 



0.36±g:?J 

17+0.05 
u - ± '-0.02 
01 o+0.004 
u.uio_ 002 

r ccr+0.06 
°- oo -0.05 

5.4 ±0.1 



0.5111°; 
0.5261°,; 
1.0261°,; 
0.1931° 
0.01451° 
5.5651°,; 



047 

().->() 

033 

030 

1)7') 

03<> 

I)))') 

008 

0007 

0007 

012 

()()(> 



This work 
This work 
This work 
This work 
This work 
This work 
KOW 



Note: All errors quoted are 68 per cent confidence. 
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4.2 Model optimization 

We used three types of optimization schemes to fit the light 
curves, a routine adapted from Bevington (1969) based on 
the Levenberg-Marquardt scheme, the "grid search" routine 
adapted from Bevington (1969), and the "differential correc- 
tions" program supplied with the W-D code. In our experi- 
ence a "brute force" procedure of searching parameter space 
is needed to reliably find the global \ 2 minimum, rather 
than just a local minimum. We began by using several dif- 
ferent initial guesses of parameter values for the Levenberg- 
Marquardt optimization routine, including the parameters 
found by KOW, and found a set of parameter values that 
gave a relatively small value of \ 2 ■ At this point it became 
apparent that the overall light curve is insensitive to the 
value of the polar temperature of the white dwarf-changing 
Tp ie(WD) only results in a change of the depth of the fea- 
ture near phase 0.5. We therefore fixed T po ie(WD) at its 
optimal value of 24,300 K, leaving four free parameters. 
Next, we defined a grid of values in the inclination-mass 
ratio plane covering a range of 75 < i < 90° in steps of 0.1° 
and 0.30 < Q < 2.00 in steps of 0.1. We started at the point 
in the i — Q plane corresponding to the provisional solution 
found above (i — 82.6° and Q = 0.76) and the routine was 
run at each point on the grid. For these fits, the values of 
i and Q were fixed at the values corresponding to the grid 
location, leaving only the two surface potentials as free pa- 
rameters. The grid search routine was then used to optimize 
the Sis. Figure 5 shows a contour map of the \ 2 values. We 
found the minimum \ 

value of Xmin = 94.84 for 100 data 
points occurred at the point i = 89.0° and Q — 1.15, rela- 
tively far from our provisional fit found above. The standard 
deviation of the residuals (i.e. the data minus the model fit) 
is 0.002 mag, which is comparable to the statistical errors of 
the points in the binned light curve. 

Our values of i = 89.0° and Q = 1.15 are signifi- 
cantly different from the values of i = 78.05 ± 0.50° and 
Q = 0.87 ± 0.15 found by KOW. We believe our present 
results are more reliable than those of KOW for two rea- 
sons. First, we have « 36.6 hours of CCD observations from 
a 2.1m telescope (of which ~ 8.2 hours was spent reading 
out the CCD), compared to 7.9 hours of observations ob- 
tained by KOW using a high-speed photometer on a 2.1m 
telescope. Hence our statistics will be better. Second, our 
CCD data are more reliable since we were able to identify 
and remove the extra 7.8 hour modulation. KOW were un- 
aware of the extra 7.8 hour modulation, and as a result their 
folded U and B light curves are likely to be biased and not 
accurately representative of the true ellipsoidal light curves. 



4.3 Computation of parameter uncertainties 

We used a Monte Carlo procedure to compute confidence 
limits on the fitted and derived parameters. At each point 
in the i-Q plane shown in Figure 5 we have fitted values 
of f2 a dB and fiwD and a \ 2 of the fit. (The fitted values 
of the Sis are basically independent of the inclination and 
hence are only a function of the mass ration Q-larger val- 
ues of Q result in larger values for both S7 s< jb and £>wd.) 
One can define a region in the plane corresponding to a 
certain confidence limit based on the change in the \ 2 val- 
ues, i.e. the 68 per cent confidence region is defined by 




Mass (M Q ) 

Figure 6. Top: The marginal distribution of the values for the 
total mass as generated using the Monte Carlo procedure. The 
dashed-dotted lines indicate the region that contains 68 per cent 
of the area and the dotted lines indicate the region that contains 
90 per cent of the area. Bottom: Similar to the top, except the 
marginal distribution of the values for the sdB mass is displayed. 



X 2 = Xmin + 2.34. We divided up the plane into 50 regions 
with region 1 corresponding to confidence limits between 
and 2 per cent (Xmin < X 2 < Xmin + 0.041), region 2 cor- 
responding to confidence limits between 2 and 4 per cent 
(Xmin + 0.041 < x 2 < Xmin + 0.081), and so on. The con- 
fidence limits defining the last region were 98 per cent and 
99.99 per cent, or Xmin+7.82 < X 2 < Xmin+18-43. Then 2000 
random i,Q pairs were selected from each region and the SI 
values were determined using a two-dimensional cubic spline 
interpolation routine (Press et al. 1992). Finally, for each i,Q 
pair, 10 values each of the mass function and orbital period 
were drawn randomly from appropriate Gaussian distribu- 
tions and the system properties (i.e. component masses, or- 
bital separation, etc.) were computed. Confidence limits on 
the parameters were computed from the marginal distribu- 
tions of each parameter. Figure 6 shows the frequency distri- 
butions for the total mass and for the sdB mass. The mode of 
the sdB mass distribution is at M s dB = 0.38 Mq, 68 per cent 
of the values are in the range 0.22 M < M sdB < 0.69 M© 
(the lines denoting the confidence limits intersect the his- 
togram at the same N- values), and 90 per cent of the values 
are in the range 0.15 M < M sdB < 1.30 Mq. The corre- 
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sponding values for the total mass are M tota i = 0.84 Mq for 
the mode, 0.53 Mq < M to tal < 1.34 M© for the 68 per cent 
confidence interval, and 0.39 Mq < M to tal < 2.22 M for 
the 90 per cent confidence interval. Table 1 gives the sys- 
tem parameters for KPD 0422+5421 and their 68 per cent 
confidence regions. 

4.4 Additional astrophysical constraints 

The 68 per cent and 90 per cent confidence intervals on the 
masses are somewhat large, which is not surprising consid- 
ering the relatively large range of mass ratios allowed by the 
W-D fits. There are potentially two separate astrophysical 
constraints we can use to limit the parameter ranges, which 
we describe below. 

4-4-1 The white dwarf mass-radius relation 

It is highly likely that the unseen star is a white dwarf, and 
as such has a maximum allowed mass and a well-defined 
mass-radius relation. However, the W-D code makes no as- 
sumption as to what the nature of the unseen star is. The 
mass and radius of the unseen star are adjusted as neces- 
sary when fitting the light curve. In our particular situation, 
Qwd is positively correlated with Q. Since the mass of the 
white dwarf is inversely correlated with Q and the radius of 
the white dwarf is inversely correlated with f2wo, we find 
that the mass of the white dwarf in the W-D fits is pos- 
itively correlated with its radius (see the bottom panel of 
Figure 7). The actual mass-radius relation for white dwarfs 
is in the opposite sense: the radius of a white dwarf get 
smaller as its mass gets larger. The form of the mass-radius 
relation is model dependent, and we used the evolutionary 
tracks of Althaus & Benvenuto (1997, 1998) and Benvenuto 

6 Althaus (1999). These authors have computed evolution- 
ary tracks for both CO and He white dwarfs using a range 
of hydrogen envelope masses and for a range of metallici- 
ties. For each track, we plot in the lower panel of Figure 

7 the radius of the white dwarf as a function of its mass 
when T cB » 35, 000 K and when T c ff ~ 15, 000 K. Some of 
the low-mass He models do not get as hot at 35,000 K. In 
these cases we plot the radius corresponding to the hottest 
available temperature. The intersection of the theoretical 
mass-radius band and the W-D model confidence limits is 
relatively small. The lower envelope of the theoretical mass- 
radius band intersects the 95 per cent confidence contour at 
Mwd = 0.42 Mq and the upper envelope intersects the same 
contour at M W d = 0.64 Mq. W-D fits with M W d & 0.4 M 
or with Mwd ^ 0.7 Mq can be firmly ruled out since the fit- 
ted radii are well outside the theoretical mass-radius band. 
If we use a "filter" in our Monte Carlo simulation and reject 
the computed values when the white dwarf mass and radius 
fall outside the theoretical mass-radius band, we find the 
following: Af sdB = 0.51lt°;°so M , M W d = 0.526±g'°^ Mq, 
and M^ai = l-026jlo'ogg Mq, 68 per cent confidence, or 
M s dB = O.51li°;°^M 0J Mwd = 0.5261°;°^ Mq, and 
Mtotai = 1.026±2; 2l M ©> 90 P er cent confidence. The mass 
of the sdB star is quite consistent with the "canonical" ex- 
tended horizontal branch mass of 0.5 Mq (Caloi 1989; Dor- 
man, Rood, & O'Connell 1993; Saffer et al. 1994). Further- 
more, the total mass of the system appears to be well be- 
low 1.4 Mq. We list in Table 1 the values of the derived 




0.2 0.4 0.6 0.8 



M WD ( M ©) 

Figure 7. Bottom: A contour plot of the two-dimensional fre- 
quency distribution of the radius of the white dwarf vs. its mass. 
The contour levels indicate the regions that contain 68 per cent 
of the volume and 95 per cent of the volume. The star marks the 
mode. We also show the theoretical mass-radius relation, taken 
from the models of Althaus & Benvenuto (1997, 1998) and Ben- 
venuto & Althaus (1999). The filled circles are for CO white 
dwarfs with 15, 000 K < T e g < 30, 000 K and with various hy- 
drogen envelope masses and the filled triangles are for He white 
dwarfs. The vertical dashed-dotted lines indicate roughly where 
the upper and lower envelopes intersect the 95 per cent confidence 
contour. Top: A contour plot of the two-dimensional frequency 
distribution of the mass of the sdB star vs. the mass of the white 
dwarf. The 68 per cent and 95 per cent confidence regions are dis- 
played. The diagonal solid line indicates Mwd + M s( jb = 1.4 Mq. 

parameters (with 68 per cent confidence limits) when the 
mass-radius constraint is imposed. 

It is quite clear that imposing the white dwarf mass- 
radius constraint greatly reduces the uncertainties on the 
derived parameters. Although the form of the mass-radius 
relation is model dependent, the constraint is fairly robust 
since the theoretical mass-radius band is nearly perpendicu- 
lar to the ridge-line of points derived from the W-D fits and 
Monte Carlo simulation. In addition, we have used models 
with a wide range of hydrogen envelope masses and metallic- 
ities, and we have also used a generous range of temperatures 
(15,000-35,000 K). 

As an aside, we can use the Althaus & Benvenuto mod- 
els to estimate the cooling age of the white dwarf. Each 
Althaus & Benvenuto model gives the mass, radius, and 
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temperature of the white dwarf as a function of the cool- 
ing age. For models where 0.50 Mq < M W o < 0.56 Mq, 
0.0138 Mq < Rwd < 0.0152 Mq, and 20,000 K < T cff < 
25, 200 K (i.e. the 68 per cent confidence ranges), we find a 
range of cooling ages between 11.3 and 50.1 million years. 
The fact that the white dwarf appears to be relatively young 
may put constraints on how the present-day binary formed. 
For example, if the white dwarf formed first, then the second 
star would have had 50 million years additional time to 
evolve into the sdB star, assuming the white dwarf has not 
somehow been heated since its formation. 

4-4-2 The radius of the sdB star 

For the second method to potentially reduce the parame- 
ter space found by the W-D fits we can exploit the fact 
that the mass of the sdB star is strongly correlated with 
its radius (S7 st jB is positively correlated with Q). Hence, the 
observed rotational velocity of the sdB star will be corre- 
lated with its radius. One must assume a rotational period 
for the sdB star, however. The usual assumption for close 
binary stars is that P ro t = Porb, i.e. synchronous rotation, 
although this may not necessarily be the case (see below). 
If the rotation is synchronous, then Vrot sin i = 80 km s -1 
when M s< jb = 0.2Mq and Vrot sin i = 125 km s -1 when 
M s dB = 0.7 Mq, where the values are computed without us- 
ing the white dwarf mass-radius constraint. Thus a measure- 
ment of V TO t sin i to better than ~ 10 km s _1 would allow one 
to reduce the available parameter space for the mass of the 
sdB star. The surface gravity computed from the W-D fits is 
also correlated with the computed mass of the sdB star. We 
find log g = 5.46 (cgs) when Af s( jB = 0.2 Mq and log g — 5.59 
(cgs) when M sa _B = 0.7 Mq, again where the values are com- 
puted without using the white dwarf mass-radius constraint. 
A spectroscopic measurement of log<? is model dependent 
and requires high quality data (one slight advantage of us- 
ing log g to constrain the range of sdB masses is that one 
does not need to assume a rotational period for the sdB 
star). Unfortunately, the spectroscopic measurement given 
by KOW of log<? = 5.4 ±0.1 does not provide a definitive 
answer. The 2a range covers nearly all of the logg range of 
the W-D fits and Monte Carlo simulation. 

4.5 A refined orbital ephemeris 

The more precisely defined eclipse profiles also allow us to 
determine the orbital phase much more accurately. The dif- 
ferential corrections routine of the W-D code can be used 
to determine the optimal phase (relative to some assumed 
time) and its probable error. The heliocentric time of the 
superior conjunction of the sdB star is given in Table 1. 
The formal la error on this time is 4.9 seconds, compared 
with the error of 43 seconds on that measurement given 
m KOW0. Our determination of the phase combined with 
KOW's determination leads to the improved measurement 
of the orbital period given in Table 1. 

The accuracy with which we can determine the orbital 
phase allows for the possibility of measuring the change in 

t Note that the error on KOW's value of To (photo) should be 
0.0005 days, not 0.005 days as given. 



the binary period. A possible cause of a change in the orbital 
period would be the loss of orbital angular momentum via 
the radiation of gravitational waves. Ritter (1986) gives the 
period derivative due to gravitational wave radiation as 

Pgr = -2.18xl0" 14 M sdB MwD(A/ sd B+A/wD)" 1/3 P" 5/3 (2) 

where the masses are in solar masses and the period is in 
days. Using the masses computed using the white dwarf 
mass-radius constraint and P given in Table 1 we find 
Pgr = -(3.lt°o.i) x 10" 13 . The phase difference A0 O _ C 
which accumulates over N orbital cycles between the ob- 
served orbital phase and the expected orbital phase based 
on a constant period is A</> D _ C « 0.5PN 2 (Ritter 1986). In 
ten years (40,502 orbital cycles) the accumulated shift will 
be « 2 seconds. The formal la error of our present phase 
determination is 4.9 seconds. One would need more than 
twenty years of phase measurements to begin to detect a 
marginal phase shift, assuming future measurements of the 
phase have similar uncertainties as ours. The phase might 
be more precisely determined if the ingress and egress por- 
tions of the eclipse profiles were better resolved, although 
this would be a somewhat daunting task as the white dwarf 
takes only « 45 seconds to move a distance equal to its di- 
ameter. Nevertheless, a measurement of Pgr would be valu- 
able since it would provide an independent constraint on 
the total mass of the system. This assumes that there is no 
other mechanism acting to remove orbital angular momen- 
tum. Note that evolutionary changes of radius of the sdB 
star on a 10 s year time-scale might lead to a P term of 
similar size to Pgr, if a small portion of the total angular 
momentum of the binary is stored in synchronized rotation 
of the sdB star. 



5 POSSIBLE SOURCES OF THE 7.8 HOUR 
MODULATION 

What is the physical cause of the 7.8 hour modulation which 
we assert is real? The white dwarf is too faint to account for 
the 7.8 hour modulation since the depth of the white dwarf 
eclipse is 0.006 mag whereas the range of the 7.8 hour mod- 
ulation is » 0.02 mag from minimum to maximum. Hence 
the extra (non-ellipsoidal) light must be from the sdB star. 
We discuss a few possibilities below. 

5.1 Star spots? 

One possibility is that the 7.8 hour modulation is due to a 
spot on the sdB star that rotates in and out of view. In this 
case, the rotational velocity of the sdB star would be sub- 
stantially slower than the velocity if it were in synchronous 
rotation: Vrot ~ 28 km s _1 compared with V ro t = 100 km 
s _1 for synchronous rotation. It should be straightforward 
to obtain a high resolution spectrum of KPD 0422+5421 
and measure the rotational broadening to see whether the 
measured rotational velocity is consistent with synchronous 
rotation (the resolution of KOW's spectrum was not suffi- 
cient to place any meaningful constraints on V I0 t)- 

The W-D fits we have computed all assume synchronous 
rotation. The W-D code allows one to compute models when 
P-ot 7^ Porb- We computed a series of W-D models with 
the parameter set to 2.16/7.8=0.277. The best-fitting 
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model had i = 89.0°, Q = 1.19, and x 2 = 94.83. The 
X 2 contours in the i-Q plane looked similar to the contour 
plot shown in Figure 5. Hence our results are insensitive to 
changes in the rotational period of the sdB star. We also 
experimented with "spot" models in the following way: We 
added a spot to the sdB star and computed several orbital 
cycles in an attempt to mimic the time series displayed in 
Figure 1. Each model time series was detrended (using an 
unspotted model) to isolate the underlying spot modulation. 
We found that we could not get a modulation with a form 
matching the shape of the modulation shown in Figure 3 
using a single spot. Several spots of various sizes placed at 
strategic locations on the star may be needed to produce the 
"sawtooth" pattern that we observe. Until a definitive mea- 
surement of V TO t sini is made that demonstrates P Iot ~ 7.8 
hours (i.e. Kotsini « 28 km s _1 rather than » 100 km s _1 
expected for synchronous rotation), a more thorough explo- 
ration of spot models does not seem warranted at this time. 

5.2 Pulsation? 

Another possibility is that the modulation could be due to 
a pulsational instability in the sdB star. Recently, a new 
subclass of sdB stars known as the EC 14026 stars has been 
identified (Kilkenny et al. 1997; Koen et al. 1997; Stobie et 
al. 1997; O'Donoghue et al. 1997). The EC 14026 stars show 
modulations with amplitudes of ^ 0.01 mag and periods 
of 400 seconds. In a few cases, these modulations have 
been identified with low-order radial and non-radial p and / 
pulsation modes (Billeres et al. 1997, 1998). Many, but not 
all, of the EC 14026 stars have detectable main sequence 
companions (e.g. Kilkenny et al. 1998 and cited references) . 
One system, PC 1336-018, is in a HW Vir type binary with 
an orbital period of 0.10 days (Kilkenny et al. 1998). It is 
not clear whether the other EC 14026 binaries are astro- 
physically "close" (periods less than a few days) or "wide" 
(separations of greater than a few AU). 

KPD 0422+5421, with T cB « 25, 000 K is much cooler 
than the EC 14026 stars (30, 000 K <> T cS <, 35, 000 K, 
Billeres et al. 1997, 1998). Although the amplitudes of the 
modes observed in the EC 14026 stars are similar to what we 
observe in KPD 0422+5421, the periods are much shorter. 
Also, the pulsations observed in the EC 14026 stars have 
waveforms that are more sinusoidal, rather than the saw- 
tooth pattern we observed. Thus if KPD 0422+5421 is in- 
deed a pulsator, it probably would be in a different class 
than the EC 14026 stars. 

Whatever its cause, the 7.8 hour modulation we observe 
in KPD 0422+5421 appears to be unique among sdB stars. 
However, we should note an important point. It is often 
the case that in reducing high-speed photometry, long term 
trends (e.g. timescales of hours) are typically removed us- 
ing polynomial prewhitening or other techniques (C. Koen, 
private communication) . Thus we encourage observers to re- 
examine their existing data and to acquire additional data. If 
the 7.8 hour modulation is a real feature of KPD 0422+5421, 
then our observation should be repeatable if given enough 
care. If the 7.8 hour modulation is not real, or not strictly 
periodic, than that too will become evident with additional 
data. 



6 SUMMARY AND SUGGESTED FUTURE 
OBSERVATIONS 

We have presented additional and extensive CCD observa- 
tions of KPD 0422+5421. We have discovered that the light 
curve of this binary star contains two distinct signals, the 
2.16 hour ellipsoidal modulation and a non-sinusoidal mod- 
ulation at 7.8 hours. This extra modulation appears to be 
intrinsic to KPD 0422+5421 and is not an artifact of the 
observation or analysis. We removed this 7.8 hour modu- 
lation from the light curve and folded the detrended curve 
on the orbital period to obtain the underlying ellipsoidal 
light curve. This folded light curve now clearly shows the 
signature of the transit of the white dwarf across the face 
of the sdB star and the signature of the total eclipse of the 
white dwarf (although the latter feature is somewhat noisier 
than the former feature). We modelled the ellipsoidal light 
curve with the W-D code and derived values of the inclina- 
tion, the mass ratio, and the fl potentials, and used a Monte 
Carlo code to compute the uncertainties in the parameters. 
The marginal distributions of the component masses as com- 
puted with the Monte Carlo code have large "tails", hence 
the parameters have relatively large uncertainties. If we ap- 
ply an additional constraint and require the white dwarf's 
radius and mass to be consistent with the theoretical mass- 
radius relation, we find the derived system parameters have 
much smaller uncertainties. Using the white dwarf mass- 
radius constraint we find that the mass of the sdB star is 
somewhat lower than the value determined KOW and is 
much closer to the "canonical" mass for extended horizontal 
branch stars. The total mass of the system appears to be 
less than 1.4 Mq with high confidence (99.99 per cent). We 
also find that the white dwarf is much hotter (and hence 
much younger) than found by KOW. The resolved eclipse 
profiles allow us to determine the phase accurately, which 
in turn may allow for the measurement of the orbital period 
derivative after a reasonable amount of time. 

Our basic picture of the KPD 0422+5421 binary is se- 
cure, i.e. it contains an sdB star with a mass close to 0.5 Mq 
paired with a nearly equally massive white dwarf. Addi- 
tional spectroscopic observations should be done to improve 
the determination of the velocity amplitude of the sdB star, 
to determine its rotational velocity, and to refine the mea- 
surements of its effective temperature, surface gravity, and 
metallicity (a good measurement of the He abundance is 
needed for a precise determination of logg). An improved 
if-velocity for the sdB star would result in a more precisely 
determined mass function, which in turn would result in 
better determined masses. Accurate values of Kotsini and 
log g could be used to reduce the uncertainty in the mass 
of the sdB star, independent of the white dwarf mass-radius 
relation. Further photometric observations should be done 
in several colours to verify that the 7.8 hour modulation is 
coherent and to see if its character depends on the bandpass, 
to better sample the ingress and egress phases of the eclipse 
profiles, and to provide additional phase determinations so 
that the change in the period can be tracked over time. 
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